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ABSTRACT 

Finding the specific nucleotides that underlie adaptive variation is a major goal in evolutionary 
biology, but polygenic traits pose a challenge because the complex genotype-phenotype relationship can 
obscure the effects of individual alleles. However, natural selection working in large wild populations can 
shift allele frequencies and indicate functional regions of the genome. Previously, we showed that the two 
most common alleles of a complex amino acid insertion-deletion polymorphism in the Drosophila insulin 
receptor show independent, parallel dines in frequency across the North American and Australian 
continents. Here, we report that the dine is stable over at least a five-year period and that the 
polymorphism also demonstrates temporal shifts in allele frequency concurrent with seasonal change. We 
tested the alleles for effects on levels of insulin signaling, fecundity, development time, body size, stress 
tolerance, and lifespan. We find that the alleles are associated with predictable differences in these traits, 
consistent with patterns of Drosophila life history variation across geography that likely reflect adaptation 
to the heterogeneous climatic environment. These results implicate insulin signaling as a major mediator 
of life history adaptation in Drosophila, and suggest that life history tradeoffs can be explained by 
extensive pleiotropy at a single locus. 
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INTRODUCTION 

Many organisms display intraspecific variation in life history traits, including differences in 
reproductive timing and allocation, development, body size, stress tolerance, and lifespan. Alternative life 
history strategies have long been hypothesized to represent adaptive responses to variable selection 
(Stearns 1992), but examples with known causal nucleotides are rare (Flatt & Heyland 2011). However, 
sequence polymorphisms have been identified recently for flowering time in Arabidopsis thaliana 
(Mendez-Vigo et al. 201 3) and adaptation to freshwater environments by marine stickleback fish, which 
includes changes in life history as well as in physiology and morphology (Jones et al. 2012). One reason 
that identification of molecular targets of selection is valuable is because it can illuminate mechanisms of 
multi-trait correlation. For example, the gene Catsup in Drosophila melanogaster affects longevity, 
locomotor behavior, and bristle number, but individual polymorphisms within the locus act on the traits 
independently and do not show pleiotropic effects (Carbone et al. 2006). Likewise, coat color in deer mice 
is a multi-phenotype trait associated with the Agouti locus, but mutations within the gene appear to have 
been targeted independently by selection with minimal pleiotropy (Linnen et al. 2013). These examples 
suggest that fitness-related traits, even those affected by the same gene, may be functionally 
independent such that natural populations might achieve any genetic and phenotypic combination. 
Alternatively, life history traits routinely exhibit strong genetic correlations, including negative associations 
between lifespan and reproduction and positive associations between lifespan and stress resistance 
(Reznick 1 985; Stearns 1 991 ; Partridge et al. 2005; Vermeulen & Loeschcke 2006; Harshman & Zera 
2007; Toivonen & Partridge 2009). These associations mirror pleiotropic effects of lab-derived mutations 
and can be difficult to break in artificial selection experiments (Leroi et al. 2005; Patridge et al. 1999; 
Anderson et al. 201 1 ; but see Khazaeli & Curtsinger 201 0; Khazaeli & Curtsinger 201 3). Consequently, 
pleiotropy, and specifically antagonistic pleiotropy, in which a genetic element encodes both positive- and 
negative-fitness phenotypes, remains an important component in the discussion of how life histories 
evolve (Williams 1 957; Flatt & Heyland 201 1 ). 

D. melanogaster are distributed across environments that range from temperate to tropical, and 
exhibit variation in life history that appears adaptive. At high latitudes, populations exhibit higher incidence 
of reproductive diapause, larger body size, higher cold stress tolerance, and longer lifespan, but lower 
fecundity, relative to low latitude populations (Capy et al. 1 993; Mitrovski & Hoffmann 2001 ; de Jong & 
Bochdanovits 2003; Schmidt et al. 2005a; Trotta et al. 2006). There is substantial genetic variance for 
these traits and pervasive genetic correlations among them, indicating that selection in the local 
environment may act on some phenotypes but drive expression of others through tradeoffs (David 1975; 
Anderson et al. 2003; de Jong & Bochdanovits 2003; Schmidt et al. 2005b; Rako et al. 2007; Schmidt & 
Paaby 2008). Tolerance to environmental stress may be especially important, and the ability to resist 
desiccation, starvation and temperature stress correlates with climatic environment (Hoffmann & 
Harshman 1 999; Hoffmann et al. 2001 ; Hoffmann et al. 2005; Hoffmann et al. 2007). This framework 
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suggests a hypothetical selection regime: at high latitude, cold winters impose seasonal stress and favor 
genotypes that confer stress tolerance and overwintering ability; correlated traits such as larger body size, 
longer lifespan, slower development and lower fecundity may evolve as co-adapted responses to the 
same selection regime or by indirect selection via pleiotropy (Paaby & Schmidt 2009). 

Few loci have been found to explain the observed genetic variance for D. melanogaster life 
history (De Luca et al. 2003; Carbone et al. 2006; Paaby & Schmidt 2008; Schmidt et al. 2008; Bergland 
et al. 2012; Remolina et al. 2012, Sgro et al. 2013). Genes in the Insulin/lnsulin-like Signaling (IIS) 
pathway are good candidates because experimental manipulation of IIS mirrors the life history tradeoffs 
observed in natural populations (Partridge & Gems 2002; Tatar et al. 2003; Giannakou & Partridge 2007; 
Gronke 2010). We previously evaluated sequence variation at the Insulin-like Receptor (In R) in wild 
populations of D. melanogaster and observed striking reciprocal dines in the northern and southern 
hemispheres for the common alleles of a complex amino acid indel polymorphism (Paaby et al. 201 0). 
The polymorphism disrupts a region of glutamine-histidine repeats in the first exon (Guirao-Rico & 
Aguade 2009). In both North America and Australia, the most common allele is at low frequency in 
tropical and subtropical populations and increases in frequency with latitude, showing highest frequency 
in temperate populations; the second most common allele exhibits the inverse dine. These patterns 
appear to be responses to similar but independent selection pressures between the continents: fly 
populations in North America and Australia were founded at different times and from different source 
populations (Bock & Parsons 1 981 ; David & Capy 1 988), and nucleotide variants on either side of the 
polymorphism show neutral patterns across geography. 

Initial tests showed that the two common alleles, which differ in length by two amino acids and are 
herein designated lnR short and InR 10 " 9 , demonstrated significant effects on fitness traits. lnR shcrt , the allele 
common at high latitudes, was associated with greater stress tolerance; InR 10 " 9 , which is common at low 
latitudes, was associated with greater fecundity (Paaby et al. 201 0). Here we show that the allele 
frequency dine in North America persists five years later and that temporal changes in allele frequency 
within a single population correlate with fluctuations in the seasonal environment, which mimic the 
temperature and resource availability differences associated with latitudinal climate. We also find that the 
InR locus in general exhibits elevated patterns of clinality and seasonality relative to the rest of the 
genome. We more thoroughly tested the functional effects of the complex polymorphism by measuring 
levels of IIS and a spectrum of life history phenotypes in lnR sh0li ar\6 InR.' 0 " 9 lines with randomized genetic 
backgrounds. We find that the alleles associate predictably with the tested traits, suggesting that the 
complex amino acid indel polymorphism at InR is an important target of contemporaneous selection in 
wild D. melanogaster populations, is highly pleiotropic, and contributes to observed tradeoffs in D. 
melanogaster life history. 
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METHODS 

Allele frequency estimates. As described in Bergland et al. (In Review), between 50-200 
isofemale lines were established from six populations spanning Florida to Maine along the North 
American east coast (latitudes 25.5°N, 30.9°N, 33.0°N, 35.5°N, 39.9°N, 44.1°N). The population from 
Linvilla Orchard in Media, PA (39.9°N, 74.4°W) was additionally sampled in the spring and fall in 2009, 
201 0, and 201 1 . Within five generations a single male from each line was pooled into population samples, 
which were sequenced on an lllumina HiSeq2000 using 100-bp paired-end reads. Reads were aligned to 
the D. melanogaster reference genome v5.39 using bwa v0.5.9-r16 (Li & Durbin 2009). The indel 
polymorphism at InR was called using UnifiedGenotyper from GATK v3.1-1 (McKenna et al. 2010), which 
identified multiple discrete polymorphisms within the complex indel haplotypes that we previously 
characterized (Paaby et al. 2010). That characterization identified many haplotypes of different lengths, 
two of which were common (herein lnR short and lnR long ); the reference genome to which the new pooled 
samples were mapped has a haplotype of intermediate length (previously referred to as InR 251 ). Thus, the 
complex haplotypes, previously identified by their lengths, were decomposed into sets of shared or 
distinct alleles at discrete polymorphic sites. 

We used our prior Sanger sequence data (see GenBank accessions GQ927177-GQ927244) to 
assign alleles of the discrete polymorphisms to the lnR short and InR 10 " 9 allele classes. Here we report those 
classifications (polymorphisms are denoted by their position in the D. melanogaster reference genome 
v5.48). Three discrete polymorphisms in the pooled sequencing data discriminate unambiguously 
between I 'nR short and InR 10 " 9 : 17405631-7, 17405651-4 and 17405784. In each case, the reference 
genome allele matches lnR short . Only one of the discrete polymorphisms, SNP 17405614, was variable 
within an allele class; the alternate and reference alleles are both associated with lnR short such that the 
association between the discrete allele identities and the larger haplotypes in which they are embedded 
(I nR short and InR 10 " 9 ) are approaching equilibrium (R 2 =0.45). SNP 17405619 segregates at low frequency 
and was invariant in our Sanger sequencing data such that the lnR short and InR 10 " 9 haplotypes we 
observed shared the same discrete allele. Both the alternate allele of SNP 17405631 and the alternate 
allele of 17405634-7 match both lnR short and InR 10 " 9 . For polymorphism 17405637, two alternate alleles 
were called; the 6-nt insertion matches InR' 0 " 9 , but the 3-nt insertion is confounded with the alternate 
allele of 17405634-7 (which is called consistently). Thus, we ignored the 3-nt insertion allele at 17405637. 

Allele frequencies of SNPs across the entire InR locus and across other sites on chromosome 3R 
were called with CRISP (Bansal 201 0) as described in Bergland et al. (In Review). To test whether the 
InR locus (3R:1 7385970-1 7455043) demonstrates elevated clinality or seasonality relative to the rest of 
the genome, we calculated the proportion of SNPs within InR that were identified to be clinal (with 
associated clinal FDR < 0.01 ) or to vary repeatedly among seasons over the course of three years (with 
associated seasonal FDR < 0.5) using estimates of clinality and seasonality reported in Bergland et al (In 
Review). We note that the high FDR for seasonal SNPs was chosen to increase power for detecting 
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fluctuations at InR outside of the /nR to " 9/s ' ,ort polymorphism, which is not included in this scan or in the 
dataset evaluated by Bergland et al. (In Review). We contrasted the proportion of clinal and seasonal 
SNPs within the InR locus to 1000 random chromosomal regions matched for length, chromosome, and 
inversion status for ln(3R)Payne. 

Genotyping. To genotype flies for the experimental assays, the identity of the complex InR indel 
was determined by measuring its amplified fragment length following polymerase chain reaction (PCR) 
with a fluorescent-tagged primer on an Applied Biosystems 3100 capillary sequencer. PCR conditions, 
including primer sequences, are described in Paaby et al. (2010). That paper reported the distribution of 
alleles across geography and referred to the alleles by their PCR fragment lengths, which we erroneously 
stated as 248 bases for the high-latitude allele and 254 bases for the low-latitude allele. The actual 
fragment lengths using the cited primers are 5 bases longer. To avoid confusion, in this paper we refer to 
these two alleles as lnR short (previously InR 248 ) and InR 1 ™ 9 (previously InR 254 ). It should be noted that 
although fragment-length genotyping does not discriminate between the SNP alleles we characterized in 
or near the complex indel polymorphism (described above), only one of those SNPs (17405614) appears 
to be independent of the lnR short and InR 10 " 9 allele classes. All evidence to date indicates that fragment- 
length genotyping unambiguously discriminates between consistent, but complex, haplotype classes (see 
Figure 1A-B). 

Fly stocks. To determine whether lnR sh0li and InR 10 " 9 alleles have different effects on phenotype, 
we generated stocks with lnR short or InR 10 " 9 on the third chromosome, and for which the X and second 
chromosomes were replaced so that they were isogenic within and across lines. The genome regions 
surrounding InR on the third chromosome were randomized, from lines originally genotyped as /nR short 
and InR' 0 " 9 , across lines. To generate the stocks, we selected lnR short and InR 10 " 9 parental strains, neither 
of which carried the ln(3R)Payne inversion, and for which the X and second chromosome were already 
replaced (for the X chromosome, stock 2475; for the second chromosome, stock 6326 from the 
Bloomington Stock Center; originally described in Paaby et al. 2010). The parental strains were crossed 
and the offspring permitted to recombine for four generations, so that genetically variable regions from the 
parental third chromosomes were distributed as randomized blocks across lines. Fourteen isogenic 
stocks were subsequently established per allele by homozygosing the third chromosome via balancer 
extraction and genotyping the InR indel polymorphism. This entire scheme was performed twice, using 
lines derived from two independent populations: Mount Sinai, NY (40.95°N, 72.84°W) and Bowdoinham, 
ME (44.01 °N latitude, 69.90°W longitude). 

I 'nR sh ° n and InR 10 " 9 alleles were evaluated as homozygotes and heterozygotes, as well as with InR 
mutant and wild-type alleles derived from a laboratory strain. The InR mutant was the hypomorphic 
lnR p5545 allele generated by P-element insertion in exon 1 (Tatar et al. 2001 ) balanced over the TM3 
chromosome (stock 1 1661 from the Bloomington Stock Center, ry 506 P{PZ}lnR p5545 /TM3, ry RK Sb 1 Ser 1 ). 
The TM3 balancer contains a wild-type InR allele (herein InR™ 3 ) that is intermediate in length between 
lnR short and InR 10 " 9 : it is three nucleotides longer than lnR short , corresponding to an additional histidine, and 
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three nucleotides shorter than InR 10 " 9 , which contains an additional glutamine. 

Flies for assays. Once the stocks were established, we used three methods to generate flies for 
the phenotype assays (derived from one or both of the independent populations). The first two tested the 
effects of lnR shori and InR 10 " 3 in a wild-type background; the third tested lnR short and InR 10 " 9 effects while 
paired with either the InR hypomorphic allele or a wild-type allele in the TM3 balancer chromosome. In 
the first method, we exerted nominal control over the background variation in the stocks by crossing 14 
lnR short and 14 InR 10 " 9 stocks in a round-robin design to generate a total of 42 lines from which we 
collected flies for phenotyping (14 each of lnR short /lnR short , lnR lon9 /lnR lon9 and lnR short /lnR lon9 ). The crosses 
were conducted in standard media vials with 7 virgin females (2-5 days old) and 4 young males (3-8 days 
old). This design produced non-independent replication within each genotypic class, and because test 
stocks were never mated to themselves, also ensured that no individuals were strictly isogenic at the third 
chromosome. In the second method, we combined stocks (minimum of five) carrying the same allele 
(I 'nR short or lnR long ) but otherwise randomized variation on the third chromosome in bottles, permitting the 
stocks to continue recombining freely. To rear flies for the assays, 40 virgin females (2-5 days old) and 20 
young males (3-8 days old) were collected and mated in fresh bottles, either within or across genotype 
class, to generate the two homozygous and the heterozygous genotypes. Unless otherwise specified, the 
phenotype assays under this method were replicated via 1 0 mating bottles for each genotypic class. In 
the third method, we crossed 12 /nR short stocks and 12 InR 10 " 9 stocks to the hypomorphic lnR p5545 fTM3 
stock (in vials as in method one). From these 24 lines, we collected progeny carrying both the InR 
hypomorphic allele and the wild-type allele in the balancer chromosome for phenotyping. 

Phenotype assays. For all assays, flies were reared and assays were performed on standard 
cornmeal-molasses media at room temperature and subject to ambient light cycles. Larval density was 
kept low to limit overcrowding. We tested lnR sh0li and InR 1 ™ 9 alleles derived from Mount Sinai in all 
assays; we tested lnR shcrt and InR 10 " 9 alleles from Bowdoinham in a subset of the assays. All lnR short and 
InR' 0 " 9 alleles tested in the presence of the InR hypomorph (lnR p5545 ) or wild-type balancer (InR™ 3 ) were 
from Mount Sinai. All statistics were performed using JMPv7 (SAS Institute, Cary, NC). Within each 
analysis, we used planned comparisons to test for significant differences between lnR short and InR 10 " 9 
genotypes. 

Lifespan and lifetime fecundity . Bottle cages were populated with 30 virgin females and 30 males, 
collected over 24 hours, from lines generated by method one from the Mount Sinai population. A total of 
42 bottles, comprising 14 of each genotypic class, were maintained by changing standard media egg- 
laying plates every day (days 1-26) or every other day (days 28 onward) until all flies were dead. Bottles 
were inverted, so that dead flies collected on the media plates. Dead flies and eggs laid were scored at 
every plate change. Lifespan was analyzed by a proportional hazards model, and cumulative fecundity 
was analyzed by ANOVA. 

Early fecundity . Virgin females and males were collected over 6 hours, separated into single-sex 
vials of 5 females or 3 males each, and aged for 3 days on standard media with topical yeast. Flies were 
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transferred in mixed-sex batches (five females and three males) to media-free vials and starved for 4 
hours, then transferred again to fresh media vials and permitted to lay eggs. The number of eggs laid was 
recorded after 12 hours and analyzed by ANOVA. Flies were generated by method two from both 
populations, and also with the InR hypomorph via method three. 

Development time . Young (3-6 days old) males and females (5 flies of each sex) were transferred 
from standard media vials to media-free vials and starved for four hours, then transferred again to fresh 
media vials and permitted to lay eggs, at densities of <30 eggs per vial, for up to four hours before being 
removed. Flies derived from Mount Sinai were generated via method one; flies from Bowdoinham via 
method two, and with the hypomorph via method three. The number of adults per vial was recorded daily 
at 8 AM, 12 PM and 8 PM for all emerging adults, and eclosion time was analyzed by ANOVA. 

Body weight, lipid content and body size . Using methods and replication schemes identical to 
those described for development time, eggs were reared and permitted to eclose as adults. After 
eclosion, adults were transferred to fresh media vials containing topical yeast and aged for 24 hours, then 
frozen at -80°C. To measure dry weight, frozen flies were dried for 24 hours at 55°C and weighed to 
0.0001 g in single-sex batches of 7 individuals on a Sartorius microbalance (1 -2 batches per original vial 
cross). The batches were transferred to 1 .5 mL tubes and lipids were extracted by adding 400 uL diethyl 
ether and gently rocking the tubes for 24 hours. Flies were dried overnight at 55°C and batches re- 
weighed; lipid content was computed by subtraction. To estimate body size, left wings were removed from 
frozen flies and affixed to an index card using clear tape. Pictures of wings from a DC 300 camera (Leica) 
on a dissecting microscope were measured for L2 and L3 vein length using landmarks 6, 12 and 13 
(Gidaszewski et al. 2009) on 3-8 males and 3-8 females from each vial cross, using tpsDigv2.12 (Rohlf 
2008). By the same method thorax width was measured for the Bowdoinham and hypomorph crosses. 
Measurements were analyzed by ANOVA. 

Chill coma recovery . Freshly eclosed flies were collected over 24 hours from flies generated by 
method one from the Mount Sinai population, and sorted into standard media vials in batches of 5 males 
and 5 females each. After 3 days, males were sorted into fresh vials and aged another day. To induce 
chill coma, vials were covered in ice and placed at 4°C for three hours, then restored to room 
temperature. Time to recovery (transition to the upright position) was recorded using a video camera and 
analyzed by ANOVA. 

Cold shock and starvation . Freshly eclosed flies were collected over 24 hours from flies 
generated by method one from the Mount Sinai population and method three with the hypomorph, sorted 
into standard media vials in batches of 5 males and 5 females, and aged for 3 days. The 42 vials 
generated by method one were replicated twice. To induce cold shock, vials were inverted and exposed 
to -20°C (air temperature) for 25 minutes. After 24 hours, individuals were scored as dead or alive. To 
measure starvation resistance, flies were transferred into media-free vials containing a cotton ball 
saturated with 2 mL of water. At 51 hours, males were scored as dead or alive; at 64 hours, females were 
scored as dead or alive. Both assays were analyzed by nominal logistic regression modeling the log odds 
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of mortality/survivorship. 

Heat shock . Freshly eclosed flies were collected over 24 hours and sorted into mixed-sex 
standard media vials of 5 males and 5 females each. Flies were generated by method two from the Mount 
Sinai population (15-20 replicate vials pergenotypic class) and by method three. Flies were aged for 4 
days and then transferred to fresh media vials. Vials were inverted, placed at 25°C, and warmed to 37°C; 
after reaching 37°C (approximately 20 minutes), flies were kept at 37°C for 25 minutes and then removed 
to room temperature to recover. After 24 hours, flies were scored as dead or alive and analyzed by 
nominal logistic regression modeling the log odds of mortality/survivorship. 

Quantitative PCR. To test whether InR alleles affect IIS, we used quantitative polymerase chain 
reaction (qPCR) to determine the relative abundance of 7 transcriptional targets of dFOXO, a central 
transcription factor within the IIS pathway (Junger et al. 2003; Puig et al. 2003; Wang et al. 2005; Casas- 
Tinto et al. 2007; Flatt et al. 2008; Vihervaara & Puig 2008; Matilla et al. 2009). Total RNA was prepared 
from flies generated by method two from flies derived from Mount Sinai. Thirty virgin females from each of 
the three genotypic classes were collected over 24 hours, aged on standard media, and snap frozen. 
RNA was extracted using RNeasy (Qiagen) and reverse transcribed using iScript cDNA synthesis kit (Bio- 
Rad), following the manufacturer's protocol. Relative abundance of transcript levels was determined using 
an ABI 7500 Fast Real-Time PCR machine and SYBR Green PCR Master Mix (Applied Biosystems) by 
the DDC T relative quantitation method. Four technical replicates were used for each sample and relative 
abundance was normalized by using GAPDH2 as an endogenous control. Other work has demonstrated 
that GAPDH2 is an appropriate control for measuring relative levels of IIS (Hwangbo et al. 2004, Flatt et 
al. 2008). Primer sequences for the qPCR reaction are: GAPDH2 forward, 

GCGGTAGAATGGGGTGAGAC; GAPDH2 reverse, TGAAGAGCGAAAACAGTAGC; 4E-BP forward, 
GAAGGTTGTCATCTCGGATC; 4E-BP reverse, ATGAAAGCCCGCTCGTAGA; l(2)efl forward, 
AGGGACGATGTGACCGTGTC; l(2)efl reverse, CGAAGCAGACGCGTTTATCC; ac76e forward, 
CAG G ATG AATG ACG CCCTTTCG G ; ac76e reverse, ATGG ACACAACACATG CCAG CAGC ; dLip4 
forward, GATAGCAATGTGCGGTTGGA; dLip4 reverse, TCATCCGTCTCCAAGGTGTG; InR T1 forward, 
CACAAGCTGGAAAGA AAGTGC; InR T1 reverse, CAAACACGTTTCGATAATATTTTTCT; InR T2 
forward, GCCTCGCACTTTGCTTATGT; InR 72 reverse, AAAAACAACGACAGCGACAA; InR T3 forward, 
TTACG CCACTG CATTCGTTC; InR T3 reverse, ATGGCCTCTCTCTCCGTCTC. 

RESULTS 

Clinal and seasonal estimates of allele frequency. To test whether InR allele frequencies 
change in association with the climatic environment, we identified multiple discrete polymorphisms within 
and near the complex indel polymorphism from pooled population samples. We identified all variants 
(SNPs or small indels) that comprise the complex indel itself and three SNPs -11,-16 and +1 30 
nucleotides from the indel that segregate at frequencies above 5% (Figure 1 A-B). Four of these discrete 
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FIGURE 1. From pooled sequencing data, we identified eight discrete polymorphisms associated with the 
complex indel in the first exon of InR (A). Each polymorphism is labeled according to position on chromosome 
3R relative to the published D. melanogaster reference genome (v5.48). The discrete polymorphisms describe 
the complex indel haplotypes we previously characterized, which include two alleles, lnR short and lnR lo "3, that 
are common at high and low latitudes respectively (B). Polymorphisms that unambiguously discriminate 
between the lnR short and lnR l0 "3 alleles are denoted with a large dot; in each case, the allele present in the 
reference genome identifies lnR short . (See Materials & Methods for relationships between the other discrete 
polymorphisms and lnR short and lnR l0 "3.) All eight of the discrete polymorphisms demonstrate some degree of 
seasonality, but four exhibit especially strong patterns in which the plotted allele increases in frequency by 
approximately 20% over the winter for three consecutive years (C). Alleles assigned to the lnR shori allele class 
(large dots) generally show higher frequency in the spring and lower frequency in the fall, suggesting that 
overwintering may impose selection pressures consistent with high latitudes. As in our earlier report (Paaby et 
al. 2010), alleles assigned to the lnR shori allele class also increase in frequency with latitude (D). The plotted 
alleles are all reference genome alleles. 
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polymorphisms show clear seasonality, with repeating fluctuations in allele frequency of approximately 
20% between spring and fall time-points over three consecutive years (Figure 1 C). The other four discrete 
polymorphisms are suggestive of seasonality: they demonstrate either comparable changes in amplitude 
but a deviation in the frequency pattern at a single time-point, or a repetitive pattern at low amplitude. The 
discrete polymorphisms also show dines in allele frequency across latitude (Figure 1 D). We were able to 
unambiguously assign allele identities for three of the discrete polymorphisms to the lnR sh0li or InR 10 " 9 
allele classes (Table S1 ), which exhibited strong allele frequency dines five years previously (Paaby et al. 
2010). In each of these cases, changes in allele frequency across latitude replicated our earlier 
observations, and seasonal fluctuations in frequency for these discrete polymorphisms support the 
hypothesis that alleles favored at high latitude are favored during the winter. 

Compared to the rest of the genome, the InR locus exhibits elevated patterns of allele frequency 
change across latitude and seasonal time. The region surrounding InR contains significantly more clinal 
polymorphisms relative to other regions of chromosome 3R matched for size and inversion status (42% 
vs. 20%, respectively; Pr(%lnR > %control) = 0.986). Relative to matched control regions, the InR locus 
also contains slightly more seasonal polymorphisms (6% vs. 4.7%; Pr(%lnR > %control) = 0.8). Previous 
examination of linkage disequilibrium across InR (Paaby et al. 201 0) indicates that the complex indel 
polymorphism is largely independent of variation elsewhere in the gene, which implies that many sites 
within InR may be targets of spatially and temporally varying selection. 

Allele effects. We evaluated the effects of InR alleles on levels of IIS and multiple phenotypes 
and found that individuals carrying the lnR short allele, which is which is prevalent in high-latitude, 
temperate environments, exhibited lower levels of IIS, better survived cold temperature stresses and 
starvation, and showed sex-specific evidence for longer lifespan. Flies carrying the InR.' 0 " 9 allele, which is 
common in low latitude, warm environments, were associated with increased signaling and exhibited 
higher rates of fecundity, larger body mass, faster development time and better survival of heat shock. 
These phenotypic differences are consistent with experimental reduction of IIS as well as life-history 
tradeoffs observed in natural populations, and suggest that a single locus can act pleiotropically on a 
suite of fitness-related traits. We replicated some of the assays by using lnR short and lnR long alleles derived 
from two independent populations (Bowdoinham and Mount Sinai); the results are reported separately. 

Levels of insulin signaling. To test whether the observed phenotypic differences between 
lnR short and InR 10 " 9 are mediated by different levels of IIS, we used qPCR to measure relative abundance 
of seven transcriptional targets of dFOXO, a central transcription factor in the IIS pathway that is 
repressed by InR activity. Decreased InR activity increases target abundance, and transcripts for five of 
these (4E-BP, l(2)efl, ac76e, dLip4 and InR 77) were highest in lnR short , intermediate in the heterozygote, 
and lowest in InR 10 " 9 ; all but 4E-BP showed statistical significance (p<0.05) between /nR s " ort and InR 10 " 9 
(Figure 2). This pattern is robust considering that the lnR shcli genotype exhibited at most twofold higher 
levels relative to InR 10 " 9 : a difference that may be biologically meaningful, and realistic for wild-type 
genotypes, but close to the limit of detection by qPCR methods. Two other targets, transcripts of InR 
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itself, show no difference in abundance level between the homozygote genotypes, but significantly higher 
abundance in the heterozygote. 



r 



FIGURE 2. Relative abundance of seven 
transcriptional targets of dFOXO, a central 
transcription factor in the IIS pathway that is 
repressed by InR activity. Increased abundance 
of these targets indicates reduced insulin 
signaling; transcript abundance for l(2)eft, 
ac76e, dLip4 and InR T1 were significantly 
greater (p<0.05) in lnR short samples compared to 
InRiong these results are consistent with the 
expectation that lnR short is associated lower IIS 
than lnR l0 "3. Transcripts were measured from 
flies derived from the Mount Sinai population. 
Error bars show 95% confidence intervals. 
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Fecundity. Surprisingly, the homozygous InR genotypes showed no significant differences in 
total fecundity: a fly carrying the InR 10 " 9 allele laid on average only 5.3 more eggs over her lifetime than a 
fly carrying lnR short (p=0.6949, Table 1 ). This was unexpected because we routinely observed 
substantially greater population numbers in our InR 10 " 9 bottle cultures, and previous results showed that 
the InR 10 " 9 genotype is more fecund (Paaby et al. 2010). Examination of the eggs laid per day, however, 
revealed that the InR.' 0 " 9 genotype lays more eggs early in life (Figure 3A). We performed another assay 
to explicitly measure early fecundity, or how many eggs young mated females laid upon first access to 
fresh food. Here, we found that the lnR long genotype is nearly six times more fecund than lnR short in the 
first 12 hours of egg-laying (p<0.0001, Table 1, Figure 3B). This result was again observed in a replicate 
assay using flies derived from the second population (p=0.0005) (Table 1, Figure 3C). 

Development time. From egg to emerging adult, flies carrying the InR 10 " 9 allele developed faster 
than flies carrying lnR short . This result was replicated in flies derived from two populations: in the first, 
males and females carrying InR' 0 " 9 developed an average of 6.0 hours faster than those carrying lnR short 
(p<0.0001 ); in the second population, development time was not significantly different among male 
genotypes, but InR 10 " 9 females emerged an average of 18.2 hours ahead of lnR short females (p=0.0029) 
(Table 1 , Figure 4). 

Body weight and size. Because reduction in IIS reduces body size (Clancy et al. 2001 ; Tatar et 
al. 2001 ), we predicted that individuals carrying the InR 10 " 9 allele would be larger than individuals carrying 
lnR short . Indeed, the average dry body weights of InR' 0 " 9 flies were heavier than lnR short flies. This was true 
for alleles tested from both populations: InR 10 " 9 females and males derived from Bowdoinham were 26.1% 
and 18.6% heavier, respectively (p<0.0001 for both), while InR 10 " 9 females and males derived from Mount 
Sinai were 1 5.5% and 1 0.9% heavier (p=0.0002 and p=0.0004) (Table 1 , Figures 5, S1 ). The higher mass 
associated with the InR.' 0 " 9 genotype is partially explained by lipid content. For both populations, the 
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A Lifetime fecundity 



B Early fecundity (Mount Sinai) 
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FIGURE 3. Average eggs laid per female over lifetime 
and in the first 12 hours after mating. Lifetime 
fecundity was measured in flies derived from the 
Mount Sinai population and was not significantly 
different between any genotypes except between 
lnRshon anc | tne heterozygote (A). However, lnR l0 "3 
females laid many more eggs in the first 12 hours 
after mating than lnR short females, in comparisons 
from both Mount Sinai (DF=1, 21; F=22.79; p=0.0001) 
and Bowdoinham (DF=1, 27; F=15.55; p=0.0005) (C). 
Error bars show 95% confidence intervals. 
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FIGURE 4. Average development time, from the time 
the egg was laid to eclosion of the adult. Flies derived 
from the Mount Sinai population showed no 
significant differences in development time by sex, 
and lnR lon s flies developed faster than lnR short flies 
(DF=1 , 3006; F=1 1 3.23; p<0.0001 ) (A). Flies derived 
from Bowdoinham showed a significant sex effect; 
InRiong females developed faster than lnR short females 
(DF=1 , 244; F=9.07; p=0.0029) (B), but male 
genotypes were not significantly different (DF=1 , 278; 
F=1.42; p=0.2341) (C). Error bars show 95% 
confidence intervals. 
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InR 0 " 9 allele was associated with greater lipid mass than lnR s ° (Figures 5, S1 ); however, only alleles 
derived from Bowdoinham demonstrated a statistically significant difference when non-lipid mass was 
included as a covariate (p=0.0066 for females, p<0.0001 for males) (Table 1 ). Wings of InR 10 " 9 flies were 
also larger than those for InR 3 * 1011 , at least in the Bowdoinham population. In that comparison, the average 
L2 wing vein length was 0.24 mm longer in lnR.' ong females (p<0.0001 ) and 0.12 mm longer in InR 10 " 9 
males (p<0.0001) (Table 1, Figure 5). In flies derived from Mount Sinai, lnR sh0li females had the larger 
wings, but only by 0.02 mm (p=0.0029), and lnR short and lnR lon9 males were not different from each other 
(p=0.5999) (Table 1 , Figure S1). 
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FIGURE 5. Average dry body weight, lipid weight, length of the L2 wing vein, and wing:thorax ratio for flies 
derived from the Bowdoinham population. lnR short flies had lower dry weight (A) and lower lipid content (B) than 
InRiong f|j es for botn females (DF=1 , 37; F=1 03.96; p<0.0001 for dry weight; DF=1 , 36; F=1 1 .09; p=0.002 for lipid 
weight) and males (DF=1, 39; F=63.17; p<0.0001 for dry weight; DF=1, 38; F=27.36; p<0.0001 for lipid weight). 
lnRshort f|j es W ere also smaller, as measured by the L2 wing vein (DF=1 , 1 1 7; F=297.32; p<0.0001 for females; 
DF=1, 123; F=35.86; p<0.0001 for males) (C), but had a smaller wing:thorax ratio than lnR lo "a flies (DF=1, 117; 
F=35.88; p<0.0001 for females; DF=1, 123; F=5.40; p=0.0218 for males) (D). Results consistent with these, but 
of weaker effect, were observed in flies derived from Mount Sinai (Figure S1). Error bars show 95% confidence 
intervals. 
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An association between large size and InR 10 " 9 , the low-latitude allele, is predictable in terms of 
insulin signaling, but unpredictable given the observation that Drosophila are larger at higher latitude 
(Huey et al. 2000; De Jong & Bochdanovits 2003). Others have shown that proportionally larger wings 
relative to thorax size are associated with better flying in Drosophila (Hoffmann et al. 2007) and that the 
wing:thorax ratio increases with latitude (Karan et al. 1998), possibly because muscles in colder 
temperatures generate less power and proportionally larger wings can compensate for this (Gilchrist & 
Huey 2004). We measured thorax width to see if bigger wing:thorax ratios were associated with the high- 
latitude lnR short allele, but they were not. lnR short flies averaged smaller wing:thorax ratios relative to InR 10 " 9 
flies (p<0.0001 for females; p=0.021 8 for males) (Table 1 , Figure 5). For all samples in which we 
measured the L2 wing vein, we also measured the L3 vein, and observed qualitatively identical results 
(data not shown). 

Stress tolerance. Reduction in IIS confers increased resistance to stress (Giannakou & Partridge 
2007), and results for three of our four stress assays accord with these findings and our hypothesis that 
the high-latitude lnR short allele mediates increased stress tolerance via reduced signaling (Figure 6). 
Females carrying the lnR shori allele recovered from chill coma on average 89 seconds faster than females 
carrying InR 10 " 9 ; lnR short males recovered 123 seconds faster (although this comparison was only 
statistically significant for males; p=0.1448 and p<0.0001, respectively) (Table 1, Figure 6). Likewise, 
lnR short females and males exhibited 3.8 and 2.4 times the odds of surviving cold shock and 1 .3 and 1 .8 
times the odds of surviving starvation, respectively, compared to InR' 0 " 9 flies. However, the lnR short allele 
did not correlate with higher tolerance of heat shock: instead, lnR short flies showed 1 5.3 times the odds of 
dying in the assay compared to InR 10 " 9 flies (Table 2, Figure 6). All four of these outcomes are consistent 
with presumed selection pressures in natural populations, as flies in high latitude, temperate 
environments, where the lnR short allele is common, likely face cold temperatures and scarce food 
availability during the winter, while flies in low latitude, semi-tropical environments, where InR' 0 " 9 is 
common, should face thermal stress from higher temperatures. However, if levels of IIS mediate these 
stress phenotypes, the heat shock outcome provides a unique example of correlation reversal between 
signaling and stress. Nevertheless, the phenotypic correlations across these traits, and the correlations 
between the phenotypes and latitude, are patterns consistent with other examples of reciprocal tolerance 
for environmentally imposed thermal stresses (Hoffmann et al. 2001 ; Hoffmann et al. 2005). 

Lifespan. Insulin signaling has been established as a major regulator of longevity in many 
animals, including Drosophila: hypomorphic alleles of InR increase lifespan relative to wild-type controls 
(Tatar et al. 2001 ), as do other manipulations that reduce IIS (Clancy et al. 2001 ; Hwangbo et al. 2004; 
Gronke et al. 2010). Consequently, we hypothesized that lnR short might increase lifespan relative to 
InR.' 0 " 9 . Flies carrying the lnR short and InR 10 " 9 alleles showed little difference in rates of aging over the 
majority of lifespan, though lnR sh0li males exhibited reduced mortality late in life (Figure 7). The lifetime 
survivorship of lnR short males was significantly longer than that of InR 10 " 9 males, with a mortality risk ratio 
of 1 .3 for InR 10 " 9 over lnR short (p=0.0095). The female mortality risk ratio for InR 10 " 9 over lnR short was only 
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FIGURE 6. Average time to recovery from chill coma and proportion of flies surviving cold shock, starvation and 
heat shock. Although both lnR short females and males recovered faster from chill coma than did lnR l0 "3 flies, the 
effect was only significant for males (DF=1, 187; F=2.14; p=0. 1448 for females; DF=1, 207; F=16.31; p<0.0001) 
(A). A greater proportion of lnR short flies survived cold shock (DF=2; x 2 =35.70; p<0.0001 for females; DF=2; 
X 2 =12.64; p=0.0018 for males) (B)and starvation (DF=2; x 2 =12.36; p=0.0021 for females; DF=2; x 2 =7.03; 
p=0.0298 for males) (C), but more lnR'°"3 flies survived heat shock (DF=2; x 2 =55.01 ; p<0.0001 for both sexes) 
(D). These flies were derived from the Mount Sinai population. Error bars show 95% confidence intervals. 
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FIGURE 7. Survivorship curves, for females (A) and males (B), from flies derived from the Mount Sinai 
population. In both sexes the heterozygote lived significantly longer than either homozygote genotype (Table 2). 
In males, the lnR shori genotype exhibited a reduced rate of aging late in life to produce a significant lifespan 
extension relative to lnR lon 3 (DF=1 ; x2=6.73; p=0.0095), but there was no difference in longevity between lnR shori 
and lnR lo "3 females (DF=1; x2=0.02; p=0.8909). 
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1 .01 , with no significant difference in longevity between these genotypes (p=0.8909). Interestingly, the 
heterozygote of both sexes lived significantly longer (p<0.0001) than either of the homozygote genotypes 
(Table 2). We hypothesize that the heterozygote longevity is a result of heterosis across the extended 
/nR-embedded locus, a consequence of the fact that the homozygotes carried identical alleles with 
haplotypes extending beyond InR. Any general increase in fitness associated with greater heterozygosity 
was probably not mediated by IIS, as the heterozygotes were associated with intermediate levels of gene 
expression for the majority of IIS targets. 

Allele effects in an InR' background. To examine how the lnR short and InR 10 " 9 alleles performed 
in a genetic background with reduced InR activity, we crossed them into a stock with a hypomorphic allele 
(lnR p554S ) over a wild-type copy of InR on the TM3 balancer chromosome (InR™ 3 ), generating the four 
possible genotypic combinations. Previous work with this hypomorphic allele showed that a single copy of 
lnR p5545 reduces InR kinase activity by approximately half, and not as much as a double mutant (Tatar et 
al. 2001). The genotypic combinations allowed us to compare single copies of lnR short to InR.' 0 " 9 , and wild- 
type flies to flies carrying a single copy of the hypomorphic InR mutant allele. 
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Consistent with the expectation that InR function increases fecundity, flies carrying the 
hypomorphic allele laid fewer eggs in the first 12 hours after mating than wild-type flies (p=0.0109). When 
paired with the hypomorphic allele, lnR long females laid on average 38% more eggs than lnR short females, 
which is consistent with the direct allele tests but was not statistically significant in this assay, and there 
was no difference in egg-laying between lnR short and InR 10 " 9 in the wild-type background (Table S2, Figure 
8). As in the direct allele tests on Bowdoinham flies, here patterns of development time were sex-specific. 
Females carrying the InR hypomorph developed 9.1 hours slower than InR wild-type flies (p<0.0001 ), but 
females carrying InR 311011 and InR 10 " 9 were not different from each other. Males with lnR sh0li developed 1 .3 
hours slower than InR' 0 " 9 males (p=0.0064), but males with or without the hypomorphic allele were not 
different from each other (Table S2, Figure 8). These patterns are consistent with the expected 
association of high InR activity and faster development, and the difference between the wild-type lnR sh0li 
and InR' 0 " 9 alleles in males is much subtler than the difference between the wild-type and mutant alleles 
in females. 

Dry weight of wild-type flies was greater than those with the InR hypomorph (p=0.0061 for 
females; p=0.0763 for males), consistent with the prediction that increased InR activity results in larger 
mass, but there were no differences between I : nR short and InR 10 " 9 (Table S2, Figure 9). However, single 
copies of lnR sh0li and InR.' 0 " 9 were weakly associated with lesser and greater lipid levels, as we saw in the 
direct allele tests (p=0.0729 for females; p=0.0624 for males), but there were no differences between InR 
hypomorphic and wild-type alleles (Table S2, Figure 9). The observations of lnR short with lower lipid level 
suggest that minor reduction in InR activity can reduce lipids even though the more general relationship 
between IIS and lipids is a negative one. This finding is in fact consistent with earlier characterization of 
InR mutant alleles, in which single copies of hypomorphic alleles were associated with reduced lipids 
(though not in our data), and double mutants with increased lipid content, relative to the wild-type (Tatar 
et al. 2001 ). Female wing sizes were not different from each other across any of the genotypic 
comparisons, but wings of InR 10 " 9 males were bigger than wings of lnR short males (p=0.01 79), as we saw 
in the direct allele tests (Table S2, Figure 10). The /nR stort and InR 1 "" 9 alleles had no effect on the 
wing:thorax ratio, but both sexes had bigger ratios in the mutant background relative to the wild-type 
background (p<0.0001 for females; p=0.0001 for males) (Table S2, Figure 10). The correlation between 
the hypomorph and larger ratio suggests that decreased InR activity may increase the wing:thorax size 
relationship, which is what we observed in the direct allele tests as well. 



DISCUSSION 

Our results indicate that the complex indel polymorphism in the first exon of InR is likely a direct 
target of natural selection in wild populations of D. melanogaster. Allele frequencies exhibit persistent 
latitudinal dines, change rapidly and cyclically over seasonal timescales, and the two common alleles are 
pleiotropically associated with multiple traits that correspond to expected pressures in distinct climatic 
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environments. The alleles are also associated with predicted up- and down-regulation of IIS, strongly 
suggesting that they mediate signaling behavior to effect differential trait expression. The fact that the 
alleles encode different amino acid sequences suggests that this occurs through changes in protein 
function. Overall our findings provide an example of life history tradeoffs influenced by a complex, 
strongly pleiotropic polymorphism. 

While the phenotypes associated with lnR short and lnR long were specific and predictive, our assays 
did not precisely test the effects of the glutamine and histidine changes that account for the length 
differences in the complex indel polymorphism. This is for two reasons. First, the lnR short and lnR long 
alleles themselves contain several other SNPs, each of which exhibits nonrandom patterns of allele 
frequency across latitude and seasonal time (Figure 1 ). Only one of these (a non-synonymous SNP at 
position 17405614) has the potential to confound our functional results, however, as the other SNPs did 
not vary between our lnR shcrt and InR 10 " 9 lines. Allele frequencies at 17405614 exhibited substantial 
seasonal and clinal patterns, though not as strongly as several of the discrete indel variants that comprise 
the complex indel polymorphism (Figure 1). Second, genome regions flanking the complex indel 
polymorphsim were not fully randomized due to the limits of experimentally imposed recombination. 
Consequently, in our test lines, the lnR short and InR 10 " 9 alleles were embedded in haplotypes that extended 
beyond the InR locus. Based on prior estimates of linkage disequilibrium, however, flanking variants 
should be unlinked to the focal indel polymorphism (Paaby et al. 201 0). Although these other variants 
may be targets of selection — a conclusion supported by our observation that the entire InR locus exhibits 
elevated clinality and seasonality relative to the rest of the genome — the fact that the indel polymorphism 
exhibits such strong changes in allele frequency across climatic environments provides complementary 
evidence for allele functionality. This accords with other reports of high recombination and rapid decay of 
linkage disequilibrium (Mackay et al. 2012), and evidence that selection precisely targets functional loci, 
including InR (Fabian et al. 2012), in natural populations of D. melanogaster. 

Two of the discrete polymorphisms associated with the complex indel polymorphism, 17405631 
and 17405634-7, do not discriminate between lnR short and lnR long but rather discriminate between the 
published reference genome genotype (which is of intermediate amino acid length) and I 'nR shcrt or InR' 0 " 9 . 
Both of these discrete polymorphisms show strong seasonal patterns of allele frequency, which suggests 
that the intermediate allele may be a target of selection as well. The discrete 17405631 and 17405634-7 
alleles segregate at substantial frequencies in all populations, and in both North America and Australia 
the intermediate-length haplotype is the third most common allele class (Paaby et al. 2010). In our 
previous survey, the intermediate allele showed no dine in North America and a weak dine in Australia. 
Here, the discrete 17405631 and 17405634-7 alleles again show little evidence of clinality in North 
America (Figure 1 D, Table S1 ). These observations support the conclusion that seasonal and 
geographical environments may impose different evolutionary forces, including aspects of demography 
and selection, even in the face of generally concordant responses to shared climatic selection pressure 
(Bergland etal. In Review). 
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Since flies derived from high latitudes live longer than those from low latitudes (Schmidt et al. 
2005a), and reduction in IIS extends lifespan (Clancy et al. 2001 ; Tatar et al. 2001 ; Hwangbo et al. 2004; 
Gronke et al. 2010), we predicted that the lnR short allele might be associated with longer lifespan than the 
lnR long allele. lnR short males did live significantly longer than lnR long males, but the magnitude of difference 
was small compared to other traits and there were no differences between the alleles in females. The 
comparatively weak longevity effect may reflect decoupling between lifespan and reproduction that has 
been observed elsewhere (Khazaeli & Curtsinger 2010) — yet it may also represent true functional 
variation. Laboratory-derived InR mutations have shown a dramatic effect on lifespan, but only the 
heterozygous combination of two hypomorphic alleles, of all genotypes tested, significantly extended 
lifespan (Tatar et al. 2001). This suggests that US-mediated traits may be determined by precise titration 
of IIS levels. 

In addition to longer male lifespan, the lnR short allele was also associated with better cold and 
starvation tolerance and lower fecundity in females, a classic example of antagonistic pleiotropy (Williams 
1957). The differences in fecundity were limited to early life, and represent variation in egg-laying rate, 
because total lifetime fecundity was not significantly different between the alleles. Modest differences in 
egg-laying in early life could have massive fitness consequences in the wild, particularly on patchy 
substrates with larval competition, yet the competitive advantage may be undetectable in the lab if 
fecundity is measured over more than one day. Consequently, we emphasize that estimates of 
reproduction and other quantitative fitness traits should be examined carefully, particularly in the context 
of decoupling correlated traits. We also note that although the high-latitude lnR short allele is associated 
with smaller body size, in general high-latitude flies are bigger (Huey et al. 2000; De Jong & Bochdanovits 
2003). If InR contributes to selection-mediated differences in body size, it either acts epistatically with 
other body-size loci or suffers antagonistic selection pressures across multiple fitness axes. This finding, 
like observations of the pleiotropic polymorphism at neurofibroma 1 in which the high-latitude haplotype 
is negatively associated with wing size (Lee et al. 2013), demonstrates the complexity of selection 
dynamics in natural populations. 

One hypothesis for strong genetic correlations among life history traits, including the tradeoff 
between lifespan and reproduction, is the presence of pleiotropic alleles. (Williams 1957; Flatt & Heyland 
201 1 ). If many genetic variants for life history are pleiotropic, then correlations among traits can never be 
completely dissolved. Our results implicate extensive pleiotropy, but the complexity of the InR indel 
polymorphism limits our ability to map phenotype to specific nucleotides. The two alleles we tested, 
lnR short and InR 10 " 9 , differ at four amino acids across a span of 1 6 residues (Paaby et al. 201 0). Even if in 
isolation these sites could break pleiotropy by functioning independently, over short timescales the 
distinction between pleiotropy and close linkage may be inconsequential (Paaby & Rockman 2013). 
However, this region at InR clearly experiences sufficient recombination (or mutation) over longer 
timescales to generate a hyper-variable set of alleles on which selection may act. Tests on additional InR 
alleles would help to resolve the functional roles of individual nucleotides. Moreover, our findings in 
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support of pleiotropy may not represent the majority of life history alleles segregating in natural 
populations. Other work has shown that recombination can generate genotypes with positively correlated 
effects for both lifespan and reproduction, which suggests that life history variation in nature is likely 
determined by both pleiotropic and recombining non-pleiotropic loci (Khazaeli & Curtsinger 2013). 
Sampling genotypes from nature may represent blocks of alleles that produce the commonly observed 
correlations, but which can be broken under the right circumstances. Finally, although the two major 
alleles of the InR indel polymorphism showed canonical tradeoffs in the lab, the role of these alleles in the 
natural environment is unknown. Whether and how they effect trait expression might vary across different 
environments (Fournier-Level et al. 2013), with the possibility that under some conditions, pleiotropy may 
present as conditional neutrality (Anderson et al. 201 1 ). Thus, true understanding of the role of this 
polymorphism in adaptive response, including fitness consequences and degree of pleiotropy, requires 
investigation rooted in the natural environment. 
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TABLE 1 . Statistical results for tests of fecundity, development, body size and chill coma recovery for InR genotypes in a wild-type background. 





Mount Sinai population 




Bowdoinham population 




Source 


DF 


SS 


F 


P 




Source 


DF 


SS 


F 


P 


Lifetime fecundity 


genotype 


2 


6160.87 


2.4026 


0.1038 














Early fecundity 


genotype 


2 


1469.05 


55.85 


<0.0001 




genotype 


2 


335.23 


8.02 


0.0018 


Development time 


genotype 


2 


81440.13 


270.48 


<0.0001 




genotype 


2 


2742.59 


1.24 


0.2891 




sex 


1 


271.65 


1.81 


0.1793 




sex 


1 


2449.16 


2.22 


0.1367 




genotype x sex 


2 


245.64 


0.82 


0.4424 




genotype x sex 


2 


9928.1 


4.5 


0.0115 


Dry body weight 


genotype (f) 


2 


2.57x10° 


19.51 


<0.0001 




genotype (f) 


2 


2.21x10° 


80.95 


O.0001 




genotype (m) 


2 


4.1 4x1 0" 7 


14 


<0.0001 




genotype (m) 


2 


4.28x1 0" 7 


33.52 


O.0001 


Lipid weight 


genotype 


2 


0.1242 


3.7185 


0.0282 




genotype 


2 


0.4206 


19.1552 


O.0001 




sex 


1 


0.1532 


9.1744 


0.0032 




sex 


1 


0.3181 


28.9744 


O.0001 




genotype x sex 


2 


0.0161 


0.4833 


0.6184 




genotype x sex 


2 


0.077 


3.5062 


0.035 




mass 


1 


0.022 


1.3187 


0.2539 




mass 


1 


0.1656 


15.083 


0.0002 


L2 wing vein 


genotype (f) 


2 


0.0253 


7.3131 


0.0008 




genotype (f) 


2 


0.5059 


189.483 


O.0001 




genotype (m) 


2 


0.0153 


4.1532 


0.017 




genotype (m) 


2 


0.3939 


21.5208 


O.0001 


Wing:thorax ratio 














genotype (f) 


2 


0.0627 


18.7883 


O.0001 
















genotype (m) 


2 


0.0628 


2.7051 


0.0708 


Chill coma assay 


genotype (f) 


2 


1262285.3 


5.49 


0.0048 
















genotype (m) 


2 


829073.78 


12.66 


<0.0001 
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TABLE 2. Statistical results for tests of temperature shock, starvation and lifespan for InR genotypes in a 



wild-type background. 





Source 


DF 


— 2 — 

X 


P 


Cold shock assay 


genotype 


2 


45.2572 


<0.0001 




sex 


1 


93.799 


<0.0001 




genotype x sex 


2 


3.0473 


0.2179 


Starvation assay 


genotype (f) 


2 


12.3564 


0.0021 




genotype (m) 


2 


7.0283 


0.0298 


Heat shock assay 


genotype 


2 


55.01 


<0.0001 




sex 


1 


0.05 


0.8294 




genotype x sex 


2 


0.65 


0.7228 


Lifespan 


genotype (f) 


2 


155.2047 


<0.0001 




genotype (m) 


2 


165.8458 


<0.0001 



TABLE S1. Characterization of the 8 discrete polymorphisms that are associated with, or are nearby, the 
complex insertion-deletion polymorphism and that are segregating at appreciable frequencies in the 
seasonal and clinal populations. 



Position on 3R 


Discriminates between 
/nR short and/nR ton9 ? 


Shows seasonality? 


Shows clinality? 


17405614 


no 


suggestive 


yes 


17405619 


no 


yes (low amplitude) 


yes (low magnitude) 


17405631 


no 


yes 


no 


17405631-7 


yes 


yes 


weak 


17405634-7 


no 


yes 


weak 


17405637 


no 


yes 


weak 


17405651-4 


yes 


suggestive 


yes 


17405784 


yes 


suggestive 


yes 
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TABLE S2. Statistical results for tests of InR and InR 9 in an InR' background. Allele refers to InR' 
versus InR 1 ™ 9 ; genotype refers to lnR p5545 versus InR™ 3 . 
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FIGURE S1 . Average dry body weight, lipid weight, and length of the L2 wing vein for flies derived from the 
Mount Sinai population. lnR short flies had lower dry weight than lnR'° n a flies (DF=1, 45; F=16.73; p=0.0002for 
females; DF=1, 44; F=14.99; p=0.0004 for males) (A), but were not different in lipid weight (DF=1, 44; F=0.14; 
p=0.7140 for females; DF=1, 43; F=1.41; p=0.2414 for males) (B). lnR lor >s females, but not males, had 
significantly larger wings than lnR short (DF= 1, 230; F=9.07; p=0.0029 for females; DF=1, 212; F=0.28; p=0.5999 
for males) (C). Error bars show 95% confidence intervals. 
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